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1 .  INTRODUCTION 


The  optical  turbulence  structure  parameter  is  a  primary  parameter  In  expres¬ 
sions  used  to  characterize  electromagnetic  energy  and  Image  propagation  through 
the  atmosphere.  Examples  of  these  expressions  Include  the  receiver  coherence 
diameter,  the  atmospheric  modulation  transfer  function,  the  Isoplanatlon  effec¬ 
tive  path  length,  and  the  scintillation  average  length  (Miller  and  Ricklln, 
1990) . 

A  basic  form  of  the  equation  often  used  for  computing  is  given  by  Tatarski 
(1961)  as 


where 

b  -  a  constant  -  3.2  (obtained  from  Wyngaard  (1973),  Hill  (1989), 
and  Andreas  (1988)), 

z  -  height  above  ground, 

kz 

Kh  =  -  turbulent  exchange  coefficient  for  heat, 

-  dimensionless  lapse  rate  ^  .  L  <  0  (Hansen,  1980), 

k  -  von  Karman’s  constant, 


-  energy  dissipation  rate  -  ~  xj'^  (Panofsky,  1968) 

-  ^  -  15^j  ,  L  <  0  (Hansen,  1980), 


-  real  index  of  refraction. 


L  -  Obukhov  length, 

^  -  height  derivative  of  n, 
az 

u*  -  friction  velocity. 

In  addition,  is  a  function  of  the  height  derivatives  of  potential  temperature 

and  specific  humidity,  that  is,  ^  and  As  such,  a  critical  question  would 

be  how  can  we  best  approximate  these  derivatives.  In  our  approaches,  we  assume 
that  the  environment  is  in  a  steady,  horizontally  homogenous  state  so  that 
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“  "H'  ^  Furthermore,  we  assume  that  the  derivatives  are 

expressible  in  similarity  form.  For  instance,  for  unstable  conditions,  the 
derivatives  of  potential  temperature  S,  specific  humidity  q,  and  windspeed  v  are 
written  as 


i  =  U  -  ^=1) 


3s  =  i/i  - 


dz  kz\ 


-  15.5V 

dz  kz\  l) 


(Hansen,  1980;  Businger.  1973;  Hoffert,  1979) 

where  the  integrated  forms  of  equations  (2)  through  (4)  are 


"  -  ^  I  • 


''  ‘  x{^"  (f^)  *  [  • 
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(9) 


y‘ 


X 


t 


(10) 


Zg  ~  roughness  length 


(11) 


L  = 


kg6l 


=  Obukhov  length 


(12) 


(Lumley  and  Panofsky  (1964);  Van  Boxel  et  al.,  1989) 

8*  -  temperature  scaling  length, 

q*  -  specific  humidity  scaling  length, 
u*  -  friction  velocity. 

Tvr  -  T,  (1  +  O.eiqJ. 

Tj  -  temperature  at  the  reference  height, 
q,  -  specific  humidity  at  the  reference  height, 
g  -  acceleration  due  to  gravity. 

e;  =  e*  +  o.eie^g  . 


(13) 


Furthermore,  since  these  expressions  are  based  on  the  notion  that  the  atmosphere 
is  in  steady  state  and  horizontally  homogeneous,  then  for  unstable  conditions 
appropriate  time  averages  are  on  the  order  of  20  to  30  min  (Hansen,  1991*).  This 
is  an  important  requirement;  that  is,  we  are  working  with  average  values. 
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Andreas  (1988)  expressed  ^  in  terms  of  potential  temperature  and  "absolute" 

humidity  scaling  lengths  6*  and  Q*.  Tunick  and  Rachele  (1991)  prefer  in 

terms  of  potential  temperature  and  "specific"  humidity  to  be  consistent  with 
Tatarski  (1961).  For  example  at  a  wavelength  of  0.55  /im  (Tunick  and  Rachele, 
1991) 

—  =  -7 .9  X  10'^  — -P  +  (l.97  X  ,  (14) 

dz  dz  \  T  dz) 


where 


P  -  pressure  in  millibars. 


T  -  temperature  in  degrees  kelvin  . 


There  are  several  ways  of  determining  6*  and  q*  for  equation  (12)  when 


dd 

Ti 


and 


4^  are  expressed  in  similarity  form  using  equations  (2)  and  (3).  Conceptually, 
dz 

the  simplest  approach  is  to  evaluate  equations  (5)  through  (7)  using  wind, 
temperature,  and  relative  humidity  data  measured  at  two  heights  using  sensitive, 
but  conventional  sensors  (Rachele  and  Tunick,  1991).  However,  experience  has 
shown  (Hansen,  1991*)  that  measurements  from  only  two  heights  generally  are  not 
sufficient  due  to  natural  variability  of  the  parameters  and  due  to  sensor  errors. 
Furthermore,  a  disadvantage  of  this  approach  is  that  the  measurement,  logging, 
and  processing  of  these  data  are  not  operationally  trivial.  Even  so,  this  method 
was  used  to  establish  distributions  of  H  and  L'E  for  this  study.  Therefore,  to 
clarLCy,  the  windspeed  variations  that  are  input  are  normally  distributed  with  a 
constant  value  used  for  their  variance.  The  H  and  L'E  inputs,  although  seemingly 
normal  with  respect  to  their  distributions  (see  figures  1  to  18  discussed  in 
section  6),  were  modeled  by  the  two-level  method  discussed  immediately  above. 
Later  in  section  6  we  discuss  the  standard  deviations  of  the  derived 
distributions  for  H  and  L^E. 


Another  approach,  which  is  theoretically  more  basic,  makes  use  of  turbulent  fluc¬ 
tuation  covariances  and  their  relationship  to  the  scaling  constants,  that  is. 


U 


=  - 


(14a) 


U*0*  = 


♦Frank  V.  Hansen,  1991,  personal  communication,  U.S.  Army  Atmospheric  Sciences 
Laboratory,  White  Sands  Missile  Range,  NM 
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where  u' ,  w' ,  and  q'  are  deviations  (fluctuations)  from  their  mean  values.  The 
overbar  Indicates  an  average  of  the  fluctuation  covariances. 

The  advantage  of  this  approach  is  that  it  only  requires  measurements  at  one 
height.  However,  the  disadvantage  is  that  the  sensors  and  logging  equipment 
suitable  for  making  these  measurements  are  delicate  and  sophisticated,  requiring 
highly  skilled  technicians  for  their  operation. 

Still  another  approach  for  determining  the  scaling  contents  is  based  on  the 
measurement  or  modeling  of  the  fluxes  of  sensible  heat  and  latent  heat  since 
these  fluxes  are  related  to  the  scaling  constants  8*  and  q*  as  follows: 


H  =  -CppU*e*  , 
L'E  =  -L'pU*g*  , 


where 


H  -  sensible  heat  flux, 

L'E  -  latent  heat  flxix, 

L'  -  heat  of  vaporization, 

p  -  density  of  air. 

The  friction  velocity  u*  In  equation  (15)  can  be  approximated  using  the 
similarity  relations  given  by  equations  (4),  (7),  and  (12). 

The  modeling  method  for  estimating  H  and  L^E  is  based  on  the  energy  balance 
equations  as  follows.  The  amount  of  solar  energy  reaching  the  ground  is  modeled 
considering  the  position  of  the  sun  relative  to  the  site  of  interest;  the  amount 
of  energy  scattered,  absorbed,  and  transmitted  by  atmospheric  gases  and  water 
vapor;  and  the  longwave  transmission  toward  the  ground.  Next,  one  considers 
energy  reflected  by  the  ground  surface  (albedo) ;  the  longwave  emission  from  the 
surface;  the  energy  absorbed  by  the  ground,  which,  in  turn,  is  either  stored  in 
the  ground,  used  to  evaporate  moisture  (resulting  in  latent  heat);  and  finally 
energy  that  heats  the  air  (sensible  heat).  The  advantage  of  this  approach  is 
that  it  is  heavily  model  oriented;  the  disadvantage  is  that  it  is  very  complex 
physically.  In  any  event,  the  purpose  of  this  report  is  not  to  explore  the  energy 
balance  approach  per  se.  Our  question  at  this  point  is  the  sensitivity  of  H, 

L'E,  and  V  in  estimating  .  Or,  put  a  different  way,  how  well  must  we  estimate 
H,  L'E,  and  V  in  our  model  to  provide  "acceptable"  values  of  .  The  question 
of  whether  we  can  do  this  well  enough,  using  the  energy  balance  approach,  will 
bo  addressed  in  a  separate  study. 
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In  this  study,  then,  we  write  in  terms  of  windspeed,  sensible  heat  flux,  and 
latent  heat  flux,  assuming  that  their  values  are  measured  or  modeled.  We  then 
consider  the  effect  on  if  estimates  of  H,  L'E,  and  V  are  assxuned  to  be  in 
error  or  vary  for  whatever  reason.  Furthermore,  we  assume  that  the  variations 
are  approximately  normally  distributed  as  found  in  a  separate  study  (Rachele  and 
Tunick,  1991)  using  the  first  method,  that  is,  the  two-level  method  discussed 
earlier. 

2.  PRIMARY  EQUATIONS 

We  write  (A  =  0.55  urn)  for  damp  unstable  conditions  (Tunick  and  Rachele, 

1991)  as 

=  A'e*'  +  +  CV’  , 

where 

a'  =  b  (6.241  X  -  Y-|p^  )  » 


B'  =  jb  [3.11  X  10-®]  )  ,  (18) 


d  =  b{3.8S  X  10-i®)-^jc-2/32-2/3|i  .  Y^j’N  }  .  (19) 

and  where 


b  -  constant  -  3.2, 
k  -  von  Karman's  constant  (0.4), 
z  =  height. 
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S*  -  scaling  constant  for  temperature, 
q*  -  scaling  constant  for  specific  humidity, 


We  determine  S'  and  q*  in  terms  of  sensible  and  latent  heat  fluxes  from  equation 
( 15) ;  that  is , 


H  X  10^ 

CpPU* 


(20) 


_  (L'E)  X  10^ 

^  ■  17^*  ' 

where 

H  =  sensible  heat  flux  (W  m"^)  , 

(L^E)  -  latent  heat  flux  (W  m”^)  , 

Cp  -  gas  constant  -lx  10^  (c.g.s.  units), 

^  -  density  of  moist  air  -  10"^  g  cm~®, 

u*  -  friction  velocity  (cms"^) , 

L'  -  latent  heat  of  vaporization  =  2.5  x  10^°  (c.g.s.  units). 


(21) 


Substituting  equations  (20)  and  (21)  into  equation  (16)  gives 


=  A' 


(  X  10^ 


+  B' 


HU'E)  X  loM 

c_l'p^u** 


+  C' 

V  L^p^u* 


-  ^6 


(22) 


The  friction  velocity,  u*,  is  a  function  of  H  and  the  Obukhov  length  L,  and  L  is 
functionally  related  to  V.  Hence,  u*  in  equation  (22)  can  implicitly  be  replaced 
by  V. 

From  equations  (17),  (18),  and  (19)  we  note  that  the  expressions  for  A',  B' ,  and 
C'  contain  the  Obukhov  length  L,  which,  in  turn,  is  a  function  of  H,  (L'E),  and 
u*,  that  is. 


L  = 


_ U^'r^^PyCp _ 

kg{H  -  0  .eiCpT^{L'E)  /l') 


(23) 
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where 

Tj  -  reference  level  temperature, 

T^j.  “  reference  level  virtual  temperature, 
g  -  acceleration  due  to  gravity, 

p„  -  density  of  moist  air  - 

3.  SECONDARY  EQUATIONS 
For  unstable  conditions  we  write  a  relationship  between  V  and  L  as 


P 

^  ■ 


(24) 


+  2  tan 


(25) 


where 


=  (i  - 14)' 


(26a) 


(26b) 


(27) 


(28a) 


dy  =  _ 2 _  ^  2  (28b) 

dx  k\{x  -  1)  {x  +  1)  1+x^ 
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However,  from  equation  (26a) 


From  equations 


dx  = 


ISZ 

4L= 


-3/4 


dL 


(28)  and  (29)  we  obtain 

cfv  =  -^lln/  ^  ~  \  +  2tan^3d  1  du* 

M  \x  *  if  L 


+ 


302U* 

4kL^ 


1 


.  _L_ 
+  1 


Equation  (30)  is  written  as 


dv  =  ttjdu*  +  a^dL  . 

We  use  equation  (23)  to  determine  the  differential  of  L,  that  is, 


where 


However,  since 


and 


dL  =  -^du*  *  Ik  ddy  , 
du*  ddy 


dL 

2u*Ty, 

du* 

kgSy 

dL 

u*’t^ 

_  rp 

ddt 

kgdt^ 

=  e* 

+  0.616^7* 

=  de 

'  +  o.eie^dg' 

(29) 


(30) 


(31) 


(32) 


(33a) 


(33b) 


(34) 

(35) 
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dL  = 


(36) 


icge; 


du* 


l—^ldQ*  +  O.610^dg*) 

kgQ*X 


Knowing  H  and  L'E,  and  having  the  equations  to  compute  u*,  we  can  determine  the 
differential  forms  of  d*  and  q*,  that  is, 


-dH  ^  Q*du* 

Cppu*  U* 


(37a) 


dg*  = 

L'pu*  u* 


(37b) 


4 .  CALCULATION  PROCEDURE 

For  this  calculation  procedure  we  assume  that  the  values  of  H,  L'E,  V^,  Tj.,  P^, 
and  fj.  (reference  level  relative  humidity)  are  known.  We  also  assume  that  T^, 
Py,  and  fr  are  precise.  The  only  parameters  that  vary  are  H,  L^E,  and  Vj.. 

1.  L  and  u*  are  computed  interatively  using  equations  (23)  and  (25). 

2.  Substituting  d0*  and  dq*  of  equation  (37)  into  equation  (36)  gives  dL 
in  terms  of  du*.  Solving  equations  (36)  and  (31)  gives  random  variations  in  u* 
and  L  for  random  errors  in  H  and  (L^E). 

3.  Variations  in  are  computed  using  equation  (22)  for  errors  in  H, 
L'  E  ,  and  u* . 

5.  DATA  SETS 

The  two  cases  considered  in  this  study  are  based  on  field  data  collected  at 
Davis,  California,  during  the  summer  of  1966  (Stenmark  and  Drury,  1970).  The 
Davis  field  site,  a  flat,  5-hectare  area  at  17  m  elevation  above  sea  level  is 
located  about  2  km  west  of  the  main  portion  of  the  University  of  California  at 
the  Davis  Campus,  24  km  west  of  Sacramento,  and  113  km  northeast  of  San 
Francisco.  The  data  were  taken  during  periods  when  the  surrounding  fields,  for 
the  most  part,  were  crop  covered  and  well  irrigated,  giving,  in  effect,  homoge¬ 
neous  surface  conditions  with  respect  to  temperature  and  moisture.  Advection 
effects  were  considered  to  be  negligible.  Profiles  of  wind,  temperature,  and 
specific  humidity  (moisture)  were  measured  at  nine  levels  from  25  to  600  cm.  Raw 
data  were  processed  to  give  1/2 -h  average  profiles.  Table  1  gives  the  reference 
level  (1  m)  values  for  windspeed  (V^)  ,  pressure,  (Pjr)  ,  relative  humidity  (f^),  and 
temperature  ,  as  well  as  the  derived  values  for  the  sensible  heat  flux  (H) 
and  latent  heat  flux  (L^E) ,  for  each  of  these  cases. 
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TABLE  1.  MICROMETEOROLOGICAL  DATA  FROM  DAVIS,  CALIFORNIA 


Case  1 

Case  2 

Date 

6-22-66 

6-3-66 

Time  (PST) 

1430 

1200 

Vr  (cm/s) 

447.9 

207.7 

Pj  (mbar) 

1000 

1000 

fr  (%) 

37.6 

35.6 

Tr  (K) 

294.35 

294.58 

H  (W/m^) 

66.57 

37.92 

L'E  (W/m2) 

515.43 

307.73 

6 .  RESULTS 

We  alert  the  reader  to  a  possible  pitfall  of  misinterpretation  and  application 
of  the  data  presented  in  this  report.  The  so-called  natural  variations  used  in 
this  study  were  not  determined  from  field  data;  instead  they  are  creations  on  our 
part  of  what  we  felt  were  reasonable.  In  particular  we  not  only  required  that 
the  distributions  of  the  fluctuations  be  normally  distributed,  but  we  also 
specified  the  values  of  the  variances.  (Note  however  that  the  mean  values  were 
determined  from  field  data- -see  section  5.)  For  example,  the  standard  deviation 
for  the  windspeed  distributions  was  3.33  cm/s.  For  this  study  we  did  not  try  to 
adjust  the  windspeed  variance  to  changes  in  the  magnitude  of  the  windspeed 
itself.  For  temperature  we  chose  a  standard  deviation  about  the  mean  of  1/3  of 
0.1  K.  For  relative  humidity  the  standard  deviation  used  to  generate  its  normal 
distribution  was  1/2  of  1  percent  relative  humidity.  The  standard  deviation  used 
for  the  pressure  distribution  wa.s  1/3  of  1  mbar. 

As  far  as  the  variances  for  the  derived  distributions  for  H  and  L^E  are  con¬ 
cerned,  they  (not  unexpectedly)  varied  from  case  to  case.  For  L'E  the  standard 
deviations  changed  from  8.0  to  11.3  W/m^,  for  cases  1  and  2,  respectively. 
Similarly,  the  change  in  standard  deviations  for  H  was  approximately  8.1  to  6.0 
W/m^. 

The  results  of  two  cases  are  presented  in  this  section  and  as  such  are  at  best 

representations  of  the  sensitivity  of  to  the  input  parameters,  as  specified, 
and  should  not  be  interpreted  as  real  world  solutions. 

Figure  1*  shows  the  two-level  model,  derived  distributions  of  H,  L'E,  and  U*  for 
Case  1.  Additionally,  the  normal  distributions  for  are  presented  to  show  the 
fidelity  of  the  normal  distributions  generated  and  used  as  input  parameters. 

Note  that  the  range  of  the  distribution  of  sensible  and  latent  heat  for  this  case 
is  approximately  50  W/m^.  The  distribution  for  U*  shows  a  narrow  range  and  rela¬ 
tive  insensitivity  to  variations  in  the  input  parameters  used. 


♦Figures  are  presented  at  the  end  of  the  text. 
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Figure  2  shows  the  derived  distributions  for  T*  (&*),  q* ,  and  L.  The  range  for 
these  distributions  is  approximately  a  factor  of  2.  Note  that  for  this  case, 
since  U*  is  relatively  large  and  T*  is  small,  the  mean  value  for  the  Obukhov 
scaling  length  is  quite  large,  representing  weakly  unstable  atmospheric  condi¬ 
tions  . 

Figure  3  illustrates  the  random  distribution  for  resulting  from  variations 
in  H,  L'E,  and  V^.  The  range  for  the  random  distribution  of  is  approximately 
a  factor  of  3.  Additionally  a  cumulative  distribution  for  is  shown.  It 

suggests  that  about  47  percent  of  the  time  will  have  a  value  equal  to  or  less 
than  its  mean  value. 

Figures  4  through  6  show  the  distributions  of  H,  L'E,  V^,  U*,  T*,  q*,  L,  and 
for  Case  2.  Note  here  that  mean  values  for  T*  and  q*  are  slightly  greater  in 
magnitude  than  those  from  Case  1,  and  the  mean  value  for  U*  is  small.  This 
results  in  L  values  lower  in  magnitude  for  this  case,  representing  more  moder¬ 
ately  unstable  atmospheric  conditions.  is  significantly  larger  in  magnitude 

(that  is,  on  the  order  of  10"^^)  and  its  range  is  approximately  a  factor  of  5. 

7.  SUMMARY  AND  CONCLUSIONS 

2 

In  this  study  we  wrote  C„  in  terms  of  windspeed,  sensible  heat  flux,  and  latent 

heat  flux.  Recall  that  our  goal  was  to  consider  the  effect  on  if  estimates 
of  Vr,  H,  and  L^E  were  assumed  to  be  in  error  or  contained  natural  variabilities. 
We  used  micrometeorological  data  from  Davis,  California,  to  evaluate  distribu¬ 
tions  for  H,  L'E,  and  finally  .  These  figures  were  discussed  in  section  6. 

We  found  that  by  constraining  the  variations  for  windspeed  to  a  a  -  3.33  cm/s  and 
with  a  range  of  approximately  50  W/m^  for  both  sensible  and  latent  heat  that  for 

these  two  cases,  cf,  could  be  determined  within  a  factor  of  3  or  5. 

Now  the  question  remains  as  to  what  are  acceptable  ranges  or  values  for  . 
The  answer  lies  wholly  with  their  use,  or,  that  is,  it  depends  on  the  application 

for  .  In  a  separate  study  (Rachele  and  Tunick,  1991),  we  found  that  the 
receiver  coherence  diameter,  can  vary  (in  one  case)  from  1.26  to  3.76  cm, 
causing,  in  turn,  a  significant  effect  upon  the  near-  and  far-field  slow  modu¬ 
lation  transfer  function.  Other  cases  resulted  in  similar  findings. 

In  conclusion  we  feel  confident  that  our  methodology  was  sound  and  that  our 
results  have  the  potential  to  suggest  how  well  we  must  estimate  H  and  L'E  in  our 
models . 
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Figure  1.  Case  1  -  Random  distribution  for; 

(d)  normal  distribution  for  V^. 
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2800  Powder  Mill  Road 
Adelphi,  MD  20783-1197 

Director 

Harry  Diamond  Laboratories 
ATTN:  SLCHD-ST-SP 
Dr.  Z.G.  Sztankay 
Adelphi,  MD  20783-1197 

Air  Force  Systems  Command/WER 
Andrews  AFB,  MD  20334-5000 

National  Security  Agency 
ATTN;  W21  (Dr.  Longbothum) 

9800  Savage  Road 

Ft  George  G.  Meade,  MD  20755-6000 

U.  S.  Army  Space  Technology 
and  Research  Office 
ATTN:  Brenda  Brathwaite 
5321  Riggs  Road 
Gaithersburg,  MD  20882 
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OIC-NAVSWC 

Technical  Library  (Code  E-232) 

Silver  Springs,  MD  20903-5000 

The  Environmental  Research 
Institute  of  Ml 
ATTN :  IRIA  Library 

P.O.  Box  8618 
Ann  Arbor,  MI  48107-8618 

Commander 

U,S.  Army  Research  Office 
ATTN;  DRXRO-GS  (Dr.  W.A.  Flood) 

P.O.  Box  12211 

Research  Trianagle  Park,  NC  27709 
Dr.  Jerry  Davis 

North  Carolina  State  University 
Department  of  Marine,  Earth,  6t 
Atmospheric  Sciences 
P.O.  Box  8208 
Raleigh,  NC  27650-8208 

Commander 

U.  S.  Army  CECRL 

ATTN:  CECRL- RG  (Dr.  H.  S.  Boyne) 
Hanover,  NH  03755-1290 

Commanding  Officer 
U . S .  Army  ARDEC 
ATTN:  SMCAR-IMI-I,  Bldg  59 
Dover,  NJ  07806-5000 

U.S.  Army  Communications-Electronics 
Command  Center  for  EW/RSTA 
ATTN:  AMSEL-RD-EW-SP 
Fort  Monmouth,  NJ  07703-5303 

Commander 

U.S.  Army  Communications-Electronics 
Command 

ATTN:  AMSEL-EW-D  (File  Copy) 

Fort  Monmouth,  NJ  07703-5303 

Headquarters 

U.S.  Army  Communications-Electronics 
Command 

ATTN;  AMSEL-EW-MD 

Fort  Monmouth,  NJ  07703-5303 


Commander 

U.S.  Army  Satellite  Comm  Agency 

ATTN:  DRCPM-SC-3 

Fort  Monmouth,  NJ  07703-5303 

Director 

EW/RSTA  Center 

ATTN:  AMS EL -EW- DR 

Fort  Monmouth,  NJ  07703-5303 

USACECOM 

Center  for  EW/RSTA 
ATTN:  AMSEL-RD-EW-SP 
Fort  Monmouth,  NJ  07703-5303 

6585th  TG  (AFSC) 

ATTN:  RX  (CPT  Stein) 

Holloman  AFB,  NM  88330 

Department  of  the  Air  Force 
OL/A  2nd  Weather  Squadron  (MAC) 
Holloman  AFB,  NM  88330-5000 

PL/WE 

Kirtland  AFB,  NM  87118-6008 
Director 

U.S.  Army  TRADOC  Analysis  Command 
ATTN:  ATRC-WSS-R 

White  Sands  Missile  Range,  NM  88002 
Rome  Laboratory 

ATTN:  Technical  Library  RL/DOVL 
Griffiss  AFB,  NY  13441-5700 

Department  of  the  Air  Force 
7th  Squadron 
APO,  NY  09403 

AWS 

USAREUR/AEAWX 
APO,  NY  09403-5000 

AF  Wright  Aeronautical  Laboratories 
Avionics  Laboratory 
ATTN:  AFWAL/AARI  (Dr.  V.  Chimelis) 
Wright-Patterson  AFB,  OH  45433 

Commander 

U.S.  Array  Field  Artillery  School 
ATTN:  ATSF-F-FD  (Mr.  Gullion) 

Fort  Sill,  OK  73503-5600 


Commandant 

U.S.  Amy  Field  Artillery  School 
ATTN:  ATSF-TSM-TA 
Mr.  Charles  Taylor 
Fort  Sill,  OK  73503-5600 

Commander 

Naval  Air  Development  Center 
ATTN:  A1  Salik  (Code  5012) 
Warminister,  PA  18974 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-M 

Mr.  Paul  Carlson 
Dugway,  UT  84022 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-DA-L 

Dugway,  UT  84022 

Commander 

U.S.  Army  Dugway  Proving  Ground 
ATTN:  STEDP-MT-M  (Mr.  Bowers) 
Dugway,  UT  84022-5000 

Defense  Technical  Information  Center 
ATTN:  DTIC-FDAC 
Cameron  Station 
Alexandria,  VA  22314 

Commanding  Officer 
U.S.  Army  Foreign  Science  & 
Technology  Center 
ATTN :  CM 
220  7th  Street,  NE 
Charlottesville,  VA  22901-5396 

Naval  Surface  Weapons  Center 
Code  G63 

Dahlgren,  VA  22448-5000 


Commander  and  Director 

U.S.  Army  Corps  of  Engineers 

Engineer  Topographies  Laboratory 

ATTN:  ETL-GS-LB 

Fort  Belvoir,  VA  22060 

TAC/DOWP 

Langley  AFB,  VA  23665-5524 

U.S.  Army  Topo  Engineering  Center 

ATTN:  CETEC-ZC 

Fort  Belvoir,  VA  22060-5546 

Commander 
Logistics  Center 
ATTN:  ATCL-CE 
Fort  Lee,  VA  23801-6000 

Commander 

USATRADOC 

ATTN:  ATCD-FA 

Fort  Monroe,  VA  23651-5170 

Science  and  Technology 
101  Research  Drive 
Hampton,  VA  23666-1340 

Commander 

U.S.  Army  Nuclear  &  Cml  Agency 
ATTN:  MONA-ZB  Bldg  2073 
Springfield,  VA  22150-3198 


Coiiunander 
U.S.  Army  OEC 
ATTN;  CSTE-EFS 
Park  Center  IV 
4501  Ford  Ave 

Alexandria,  VA  22302-1458 
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